Efficient method to calculate total energies of large nanoclusters 
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We present an approach to calculate total energies of nanoclusters based on first principles es- 
timates. For very large clusters the total energy can be separated into surface, edge and corner 
energies, in addition to bulk contributions. Using this separation and estimating these with direct, 
first principles calculations, together with the relevant chemical potentials, we have calculated the 
total energies of Cu and CdSe tetrahedrons containing a large number of atoms. In our work we 
consider polyhedral clusters so that in addition our work provides direct information on relaxation. 
For Cu the effects are very small and the clusters vary uniformly from very small to very large 
sizes. For CdSe there are important variations in surface and edge structures for specific sizes; nev- 
ertheless, the approach can be used to extrapolate to large non-stoichiometric clusters with polar 
surfaces. 

PACS numbers: 31.15.A-, 61.50.Ah, 68.35.Md 



I. INTRODUCTION 



Nanoscience provides an ideal platform in the search 
for novel materials with desirable and tunable properties. 
Nanocrystals (NCs) of various sizes and shapes have been 
found to exhibit a wide variety of physical and chemical 



propertiei 



I, rarely seen in bulk materials. Synthesiz- 



ing nanoparticles of a given size and shape is notoriously 
difficult and has become a key focus area due to techni- 
cally significant properties that depend on the size/shape 
of the cluster. Growth of such clusters is governed by 
both kinetics and thermodynamics^ . If a NC has a highly 
symmetric crystal structure (such as zincblende), it is 
likely that when synthesizing, the crystal will grow with 
no preferred direction of growth. On the other hand, if 
the crystal structure has a preferred axis of symmetry, as 
in the hexagonal wurtzite structure which has a unique 
polar axis, preferential growth along this axis can be ex- 
pected^ii. Energies associated with various facets (i.e., 
with different surface orientations) will also play a key 
role during growth. The total energy, which depends on 
these facet orientations, will determine the stability of a 
nanocluster having a given size and shape^s^iiiS. Calcu- 
lating these energies associated with large nanocrystals is 
a nontrivial task. Traditional "brute force" or direct first 
principles methods become quite laborious as the size of 
the cluster increases. 

First principles methods have provided reliable total 
energies for atoms and molecules; this is also true in 
solids, provided that the number of atoms in a given unit 
cell is relatively small. When this number becomes large 
(say more than a thousand), then it becomes computa- 
tionally prohibitive to carry out first principles calcula- 
tions. Naturally, it is desirable to develop alternate tech- 
niques to obtain quality total energies of systems contain- 
ing a large number of atoms. This work is focused on ob- 
taining such total energy estimates using first principles 



values calculated for a small number of well defined clus- 
ters. This approach provides crucial information needed 
to treat accurately these clusters as well as general results 
that can be used for other classes of clusters. 



II. METHODOLOGY 

A. Calculational method for clusters 

Surface energies have been extracted from first prin- 
ciples slab (total energy) calculations by several groups 
(see for example, Refs i^^d^ ). Besides that, a first prin- 
ciples wedge-shaped approach has also been proposed to 
calculate the surface energies of polar surfaces^. In the 
present work, our goal is somewhat different compared 
to the above. We will use energies obtained from first 
principles calculations on polyhedral nanoclusters to es- 
timate total energies of larger nanoclusters. The energy 
can be written as 



^tot 







where, CTq, e/3 and denote surface, edge and corner 
energies respectively, Ni is the number of atoms of type i, 
and fii is the corresponding chemical potential. The last 
term in Eq. ([T]) contains the bulk energy if the chemical 
potentials satisfy the conditions given in section II B. For 
a stoichiometric cluster the total energy is independent 
of the individual chemical potentials and is well-defined 
relative to the energy of the bulk crystal. 

In the present work, a least-squares fit for the total 
energies of several small clusters were used to estimate 
the surface, edge and corner energies (as parameters). 
If the above energies can be evaluated in a computa- 
tionally efficient way, then the total energy of a cluster 
(polyhedron), containing a (substantially) large number 
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of atoms, can be expressed algebraically as in Eq. ([T]). We 
demonstrate that using accurate density functional the- 
ory (DFT) based estimates of total energies of a few small 
clusters, it is possible to estimate the above (surface, 
edge, corner) contributions (as well as chemical poten- 
tials) and then use them (as parameters) in larger clusters 
to estimate total energies. The calculations were based 
on the local density approximation within the DFT— as 
implemented in the local orbital SIESTA codei^. Norm- 
conserving nonlocal pseudopotentials of the TrouUier- 
Martins typ&i^ were used to describe all the elements. 

For a specific shape, the cluster can be defined by one 
characteristic length iLg, so that the surface and edge 
terms can be expressed as Aa — aaP ^ Lp = bp£^ re- 
spectively, with Cq and bp being constants. Therefore, 
Eq. ID) may be written in the following form; 

i 

^^(Taac,+^ei3bi3/(. + ^c^/i^. (2) 

a p 7 

We label the energy expression on the left side of Eq. ^ 
as "termination energy" since it is the added energy due 
to the presence of surfaces, edges and corners. For ex- 
ample, for a tetrahedral structure bounded by four (111) 
surfaces, (see Fig. [T]), The characteristic size is ^Lq (^ 
being an integer and Lq the nearest neighbor distance 
along an edge), representing an edge of a triangular (111) 
facet. The tetrahedral structure includes four equivalent 
(111) surfaces, six equivalent (111) — (HI) edges, and 
four equivalent (111) — (111) — (111) corners. This choice 
enables us to work with well defined (111) surfaces, as 
well as equivalent corners, edges and surfaces. In this 
case, the termination energy in Eq. ^ turns out to be 

«f -5]7V,/iO/^' = V3aLg + 66LoJ+4cl.(3) 

i 

In later discussions, we use the above equation for both 
FCC Cu clusters and zincblende CdSe clusters. 

B. Non-stoichiometric case 

One of the examples we have chosen to study is the bi- 
nary compound CdSe. When the cluster is stoichiomet- 
ric, the number of Cd atoms will be equal to the num- 
ber of Se atoms, and the sum of the chemical potentials, 
P^Cd + /iSe = total energy per CdSc pair, can easily be 
evaluated from bulk total energy calculations. In a non- 
stoichiometric case, the number of Cd atoms is different 
from the number of Se atoms, and some energy terms, 
such as the surface energies, depend on the separate val- 
ues of ^cd or /i5e, instead of their sum. To examine this 
further, we set lower and upper bounds for the individual 
chemical potentials (pertaining to bulk CdSe) as 

<"f(Cd5e) = Mcd + MSe (4) 
= E\^l\Cd) + Etl^{Se) + /\Hf, (5) 



where Affy is the formation energy. The bulk values, 
E^^I''{Cd) and E^^I''{Se), are obtained from total energy 
calculations of pure Cd and pure Se (per atom, in their 
equilibrium structures) separately, while E^^j^{CdSe) is 
the total energy of a CdSe pair in the bulk. 

Following well known procedures'^, we can set bounds 
listed below for the chemical potentials of the individual 
species in CdSe: 

iSLt'^ (Cd) + AHf < f,cd < (Cd) (6) 
(Se) + AHf < /i5e < (Se). (7) 

The right hand side of the first inequality represents the 
fact that, ncd in the cluster must be smaller than the 
(pure Cd) bulk value E'^^l'^^Cd), since otherwise, Cd 
must phase separate. The left hand side of this inequality 
follows from the fact that for use, one can use the same 
argument and utilize Eq. ^ to obtain the following: 

< -/ise + E^:i\Se) ^ ^icd - Et:i\Cd) - AHf. (8) 

From the above arguments, it appears that we can 
only evaluate the individual chemical potentials, ^cd and 
/iSe, within the range given above. However, for a spe- 
cific family of non-stoichiometric clusters having the same 
shape, we demonstrate below that the total energy can 
be determined with no knowledge of the chemical poten- 
tials. 




FIG. 1: CdSe Zincblende (tetrahedral) structure bounded by 
four (111) surfaces terminated by Cd atoms. Small dots rep- 
resent Cd atoms, while large dots represent Se atoms. The 
characteristic size is ILq {£ being an integer, with £ — 6 shown 
here), representing an edge of a triangular (111) facet. 



We first show that the total energy of any CdSe 
(zincblende) tetrahedron, E^°l^ , bounded by similar 
(111) facets as shown in Fig. [1] can be expressed in a 
slightly different form of Eq. ([1]), i.e., 

Ef°ly{CdSe) - Ncdf^cd - Nsef^se 

= V3aLl£'^ + 6eLo£ + 4c. (9) 
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Here Lq is the nearest neighbor distance along an edge 
of the tetrahedron. 

Combining with Eq. ([4]), we obtain: 

Ef:ly{CdSe) - NseEt:l^{CdSe) 
= [Ncd - Nse)tJicd + ViaLlf + QeL^ + 4c. (10) 
However, note that 

Ncd-Nse = ]^f + ^t+l, (11) 

which results from a simple count of the atoms in a tetra- 
hedron having an edge of length iL^ {£ being an inte- 
ger and I/O the nearest neighbor distance along an edge). 
This result clearly shows that the difference in the num- 
ber of Cd and Se atoms arises from surfaces (£^), edges 
{£^) or corners {£'^), and leads to the following important 
simplification: 

Ef:ly{CdSe) - NseEti^iCdSe) 

1 3 

= [i^t^Cd + Vi<7Ll)e + i-ncd + GeLo)e + {ficd + 4c). 

(12) 

The significance of Eq. (fT^ is that, even in this non- 
stoichiometric case, it is possible to estimate the total 
energies independent of the chemical potentials. This is 
because in the above equation, the coefficients oi I 
and the constant term act as straightforward parameters 
to be estimated. Since there is no volume {L^) term on 
the right side of Eq. p^ . neither ficd nor ^se will have 
a direct effect on the final total energy to be predicted. 
Note that El^l^{CdSe) has an unambi guous value, as 
the total energy per CdSe pair in the (bulk) zincblende 
structure. Now the parameters, such as surface energy, 
edge energy, and corner energy, can be fitted using sev- 
eral, known (DFT based) total energy values from small 
clusters. Finally, in order to predict the total energies of 
large polyhedrons, the algebraic expression in Eq. ([72]) 
can be utilized as previously. 

III. RESULTS 

Test results from pure (fee) Cu as well as (zincblende) 
CdSe clusters show that this scheme is reliable to a high 
degree of accuracy. One of the significant results of the 
present study is our ability to calculate energies of nan- 
oclusters that are non-stoichiometric and that have polar 
surfaces. Furthermore, we will demonstrate that the to- 
tal energies can be evaluated independently of individual 
chemical potentials (Eq. ([5]) (O). 

A. FCC Cu 

We begin our discussion with pure fee Cu clusters, 
which are regular tetrahedral clusters having a charac- 
teristic length, L = ILq {£ being a positive integer. 




FIG. 2: The termination energy (see text for definition) con- 
tribution scaled by £^ for Cu clusters as a function of the 
inverse of the characteristic size {!/£) of the tetrahedrons. 
Y-intercept is proportional to surface energy (a), slope is pro- 
portional to edge energy (e), and curvature is proportional to 
corner energy (c). The solid squares represent the values cai- 
culated using DFT while the open triangles represent the es- 
timated values from our algebraic expression. The solid line 
is a least squares fit for the DFT points. 



Lq — 2. 545 A after geometric relaxation), bounded by 
four, equivalent (111) facets. The relevant chemical po- 
tential, ^cu , can be obtained from a fee bulk total energy 
calculation. The first principles total energies, Ef"]:^, for 
the tetrahedrons are obtained after fully relaxing a given 
cluster. 

We find it is more useful to focus on the energy term 
(-Bf°/^ — NculJ'Cu), since this represents a termination 
energy to the bulk chemical potential contribution due 
to the presence of surfaces, edges and corners. The ter- 
mination energies scaled by the square of a characteristic 
length £^, i.e., [Ef^l^ - iVc„MCu)/^^ for different sized 
clusters (^ = 2 to ^ = 8) calculated using DFT are shown 
as solid squares in Fig. [2l To estimate surface, edge 
and corner energies, total energies of a few small clusters 
{£ — 2, 3, 4, 5, 6) were used in a least squares fit according 
to Eq. The estimated cr, e and c were 0.1164 eV/A'^, 
0.2949 eV/A, and 1.0047 eV respectively. A previous 
DFT calculation has reported the (111) fee Cu surface 
energy a = 0.1213eV/ A'^ while the experimental re- 
sult for the surface energy of Cu(iii) is Q.llAAeV/A^ — . 
Hence our calculated a value is in good agreement with 
previous theoretical and experimental values. Using the 
three parameters for surface, edge and corner energies, 
we now estimate the termination energy as a function 
of cluster size, as shown in Fig. [2l The estimated total 
energies for the two larger clusters {£ = 7 and £ — 8) 
are carried out using the algebraic expression shown in 
Eq. ©. 
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Cu 




(eV) 


DFT 


algebraic 


2 


10 


-3041.43 




3 


20 


-2704.60 




4 


Q CC 
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5 


56 


-2726.96 




6 


84 


-2840.92 




7 


120 


-2981.98 


-2981.95 


8 


165 


-3139.40 


-3139.41 



TABLE I: Total energies for tetrahedral fee Cu clusters 



The energies predicted using the algebraic expression 
are in excehent agreement with the 'exact' calculations 
from first principles as evident from the results for £ = 7 
and ^ = 8 tetrahedrons in Table H] and Fig. [2l The above 
result clearly provides further support for our method of 
estimation. 



B. Tetrahedral, zincblende based CdSe clusters 
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FIG. 3: The termination energy contribution (see text) di- 
vided by for CdSe clusters as a function of the inverse of 
the characteristic size. The squares correspond to DFT en- 
ergies while the lines are least squares fits. The solid line 
corresponds to the minimal Hcd, while the dashed line cor- 
responds to the maximal ncd- The open triangles represent 
energies obtained from our algebraic calculation after fitting. 



The second system in our discussion is a non- 
stoichiometric, zincblende CdSe cluster bounded by four 
equivalent (111) facets terminated by Cd atoms. As de- 
scribed in the Methodology section and Fig. [1] the char- 
acteristic size is £Lo with Lq — 4.342A. The bulk values, 
E^^I''{Cd) and E^^I''{Se), are obtained from total energy 
calculations of pure Cd (in hep structure) and pure Se 
(in trigonal structure) separately, while Ef^l''{CdSe) is 
the total energy of a CdSe pair (in zincblende structure). 
Through explicit calculations, we obtain E^^^l^iCd) — 
-1467.17 eV]E'l^i''{Se) = -256.97 eV; E^^I'^lcdSe) = 
— 1724.84 eV. Therefore, the formation energy is calcu- 
lated to be AHf = -0.7 eV (from Eq. 

In Table [III DFT based total energies for polyhedrons 
from ^ = 2 to £ = 8 are shown. Using DFT based 
small clusters energies (for £ = 2,3,4,5,6), we have 
obtained the necessary fitting parameters involving sur- 
face energy, edge energy, and corner energy, which vary 
within the intervals 50-61 meV/A'^, 227-268 meV/A, 
and 39 — 214 meV respectively. We also see that cer- 
tain clusters undergo noticeable reconstructions, while 
others do not. Surface energy varies within the inter- 
val 50 — 61 meV/A^, showing an average of surface en- 
ergies within the reconstructed (for example £ = 4, 5) 
and non-reconstructed (for example £ = 6) structures, 
which is reasonably lower than previous DFT calcula- 
tions ( CT ~ 75 — 95 meV/A^ — ). Surface reconstructions 
have been studied before using various methods, such as 
the tight binding methodic. Here we have observed size 
dependent, surface reconstruction of the CdSe system by 
using DFT combined with the local density approxima- 
tion. 

In Fig. [3l for different Cd chemical potentials, we 
plot the scaled termination energy, {Ef^j:^ — NcdfJ-cd — 
NsetJ'Se)/^'^ by subtracting the chemical potential con- 



tributions from the total energy of the cluster as a func- 
tion of l/£. The algebraic values of Ef°j:^ obtained from 
Eq. ^ for £ = 7, 8, are indicated on the fitted curves, 
along with other values resulting from direct DFT calcu- 
lations. In general the termination energy depends on the 
choice of /itcd) as evident from Fig. [31 which is required 
to compare these surfaces with other surfaces. However, 
the predicted values for the total energy {Efl^°^) of the 
^ = 7, 8 clusters relative to other polyhedra with the 
same shape can be found using Eq. dTU]), which yields 
the same value for different values of ficd (see Table [Hi 
i.e., the algebraic value Sf"/''/^^ = -4032. 30ey for the 
£ = 7 polyhedron for all the values of ncd-) This value is 
comparable to the energy obtained from the direct DFT 
calculations, which is an advantage since it shows that 
we can obtain certain energies independent of ncd, as 
discussed earlier. 



£ (Lo) 


Ncu 


Nse 


Ef:ly/£^ (eV) 


DFT 


algebraic 


2 
3 
4 
5 
6 


10 
20 

35 
56 
84 


4 
10 
20 
35 
56 


-3920.98 
-3542.94 
-3528.42 
-3644.61 
-3821.40 




7 

8 


120 
165 


84 
120 


-4032.21 
-4263.30 


-4032.30 
-4263.30 



TABLE II: CdSe (zincblende) total energies for various cluster 
sizes. 
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IV. DISCUSSION 

For polyhedral nanoclusters of fee Cu, ealeulations 
from 10 to 165 atoms show that the energies are well 
described by this form even for small clusters. Thus the 
energies for all sizes can be determined efficiently based 
on calculations for small clusters and we propose that 
this is a useful approach for metals. 

When testing this approach for CdSe polyhedrons that 
are Cd terminated and non-stoichiometric, one encoun- 
ters energy contributions that depend on the individual 
chemical potentials, as seen from Eq. ([1]) which defines 
E^°l^ . These chemical potentials are subject to upper 
and lower bounds^'* (see inequalities ([6]), (O). However, 
it turns out that the total energies of even the non- 
stoichiometric CdSe polyhedra considered here can be 
obtained without knowing the individual chemical poten- 
tials. This is sufficient to extrapolate to large size for this 
family of clusters independent of chemical potentials and 
using only directly calculated total energies. The reason 
for the above is that the total energy of such clusters can 
be calculated from well defined bulk energies by adding 
surface, edge and corner termination terms; these termi- 
nation terms appear as mere parameters that scale with 
the dimensions of the cluster and can be estimated from 
DFT calculations of small clusters with similar topolo- 
gies. In addition, the surface, edge and corner termina- 
tion terms can be used for other clusters by including the 
chemical potentials in the way given in Eq. (jl2[) . 

For CdSe, we find an overall trend similar to that for 
Cu; however, there are deviations from a smooth curve 
for the energies as a function of size. The deviations are 
likely to be associated with reconstructions that are seen 
for certain nanocluster sizes. These reconstructions prob- 



ably originate from the changes in state occupations near 
the Fermi energy, with the largest changes apparently oc- 
curring for the edges. Also note that as the cluster size 
i — > cxD, the corner and edge contributions become small 
compared to the surface energies and thus providing a 
way of estimating the latter for large clusters. 

The total energies are described well by a least squares 
fit, carried out using the results for 5 clusters with 14 
to 140 atoms, which accurately determines the energies 
for clusters with 204 and 285 atoms. However, the di- 
vision into three separate contributions as surface, edge 
and corner energies, is not as well determined due to 
the (possible reconstruction-induced) variations between 
the different clusters. The estimated surface energy rep- 
resents an average of different reconstructed and unre- 
constructed clusters, which is reasonably lower than the 
value reported in previous DFT calculations, as pointed 
out in the results section. 



V. CONCLUSION 

In summary, we have presented an efficient method 
for calculating total energy of large nanoclusters using a 
parametrized, algebraic form with parameters fitted from 
small, first principles based, nanocluster calculations. 
The method appears to work quite well for pure met- 
als. Even for non-stoichiometric, semiconducting clus- 
ters with polar surfaces, this approach provides a way of 
estimating total energies of large clusters, important in- 
formation such as relaxation energies that are specific to 
the chosen clusters, as well as surface, edge and corner 
energies that can be used for other clusters. 
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